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LIST  OF  SYMBOLS 


received  signal  power,  watts 
transmitter  power,  watts 

transmitting  antenna  gain  in  the  target  direction 
receiving  antenna  gain  in  the  target  direction 
bistatic  radar  cross-sect  ion ,  square  meters 
transmitter  to  target  distance,  meters 
receiver  to  target  distance,  meters 
propagation  loss  over  the  transmitter  to  target  path 
propagation  loss  over  the  receiver  to  target  path 
system  loss 
wavelength,  meters 


y 

i 
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power  density  at  the  receiver,  watts  per  square  meter  ;J 

power  density  at  the  target,  watts  per  square  meter  f 

l 

unit  vector  in  the  direction  of  the  reflected  ray  3 

1 

•  r  3 

unit  vector  in  the  direction  of  the  transmitted  ray  | 

I 

distance  from  the  target  to  the  point  of  observat i on--may  3 

include  or  exclude  the  receiver  1 


unit  vector  normal  to  and  outward  from  the  prolate  spheroidal 
surface 


point  on  the  prolate  spheroidal  surface 


1.0  INTRODUCTION 


I 

I 

Early  experimental  radar  systems  were  predominantly  of  the  bistatic 
type,  where  the  transmitting  and  receiving  antennas  were  usually  separated 
by  a  distance  comparable  to  the  target  distance  [1]  -  [2].  When  the  du- 
plexer  was  developed  at  the  Naval  Research  Laboratory  in  1936,  a  means 
of  using  the  same  antenna  for  transmitting  and  receiving  made  monostatic 
(one  site)  radar  feasible.  Bistatic  radar  became  dormant  until  the  early 
1950's,  when  interest  revived  [1].  Interest  has  increased  with  the  develop 
ment  of  low  level  cruise  missiles  in  1978,  when  the  need  for  bistatic  cross 
section  results  became  apparent  for  missile  guidance  and  target  detection. 

Bistatic  radar  systems  utilize  two  separate  locations.  The  radar 
transmitter  is  found  at  one  location,  while  the  radar  receiver  is  found 
at  another  location.  in  some  configurations,  tue  radar  transmitter 
might  be  located  in  a  B~92  airplane  and  the  radar  receiver  might  be 
located  in  low  altitude  offensive  air-to-ground  missiles  such  as  the 
AGM-69A.  ,  M  other  configurations,  the  radar  transmitter  would  be  located 
t  on  the  ground  to  control  several  ground-to-air  defensive  missiles.  In 

the  first  example,  the  target  could  be  an  enemy  missile  silo  or  ground 
radar  system.  In  the  second  example,  it  would  be  either  the  attacking 
offensive  missiles  or  a  B- 52  guidance  and  control  airplane.  The  parameters 
for  detection  include  Doppler  frequency  shift  from  the  target,  signal 
frequency,  and  bistatic  radar  cross-sect  ion  created  by  the  target.  This 
study  is  limited  to  the  bistatic  cross-section  of  a  typical  target.  A 
prolate  spheroid  was  selected  because  it  resembles  an  air  cruise  missile. 
The  air  inlets  and  tail  fins  were  not  included  in  the  anal/sis  presented 
here.  Results  are  presented  in  the  form  of  bistatic  cross-sections  and 
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the  computer  program  to  generate  additional  cross-sections  for  prolate 
spheroids  with  different  major  to  minor  axis  ratios. 

2.0  BISTATIC  RADAR  CROSS-SECTIONS 

The  purpose  of  this  study  was  to  formulate  the  bistatic  radar  cross- 
section  of  a  perfectly  conducting  prolate  spheroid  when  the  prolate 
spheroid  is  large  compared  to  the  wavelength  of  the  incident  radiation. 
There  are  many  studies  [3]  “  (7]  on  the  monostatic  backscatter  from  small 
prolate  spheroids  (prolate  spheroids  whose  size  is  comparable  to  wave¬ 
length)  and  also  with  studies  of  field  distributions  at  the  boundary  of 
the  shadow  where  diffraction  effects  must  be  considered,  but  these  are  not 
considered  here.  Objects  are  almost  always  large  compared  to  wavelength 
in  detection  applications,  and  detection  is  seldom  done  near  the  shadow 
region.  In  the  region  that  we  are  concerned  with,  the  geometrical  theory 
of  optics  can  be  applied,  and  can  be  used  to  generate  some  insight  into 
the  bistatic  cross-sectional  area  as  a  function  of  the  spatial  coordinates. 


Theory : 


The  bistatic  radar  equation  can  be  written  in  the  following  manner: 


P  = 
r 


PtGt 


G>2 


4itD  2l  (r)  Ls  4-rrD  2L  (t) 
r  P  t  p 


(1) 


where  the  symbols  are  defined  in  the  prefix  of  this  report.  Eq .  (1)  is 
arranged  to  emphasize  the  motivating  analysis  in  the  derivation  of  the 
radar  equation.  The  power  density  reaching  the  target  is 

PtGt 

'Jt  ?  ’ 

4 ttD  "L  (t) 
t  P 


(2) 


the  effective  capture  area  of  the  receiving  antenna  is 


'i 

where  A“ /4 tt  is  the  universal  antenna  constant,  and  the  power  density  Ur 
at  the  receiving  site  is 


4?p  L 
r  s 


r 


(4) 


The  system  loss  is  assumed  to  occur  only  at  the  receiving  system.  Although 
it  actually  occurs  in  various  locations,  the  exact  location  or  even  its 
existence  do  not  effect  the  subsequent  analysis.  The  use  of  Eqs.  (2)  and 
(4)  can  rewrite  Eq.  (l)  as 


U 

r 


4 7t D  2l  (r) 

r  P 


(5) 


and  the  form  of  c,  suitable  for  computation  is 
o 


9 

O,  =  4ttD  L  /  1  TT— -  , 

b  r  p  (  r )  Ut 


and  to  numerical  formulation  if  Ur  can  be  found  when  is  given. 

Figure  1  shows  the  case  from  geometrical  optics  of  an  incoming  ray, 
k,  being  reflected  from  a  point  on  a  curved  surface  in  the  direction,  k'. 
A  very  important  point  arises  here  with  respect  to  the  prolate  spheroid. 
Since  the  object  is  convex,  a  one-to-one  correspondence  exists  between 
reflection  points  and  receiving  points.  As  such,  no  interference  exists 

between  refection  points,  and  U  can  be  described  simply  bv  k'  and  how 

f 

energy  density  decreases  as  a  function  of  distance  between  the  reflection 
point  and  observation  point.  In  other  words, 
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spread  (tk 1 )  v/; 

where,  spread(tk')  is  a  function  describing  the  spreading  out  of  energy 
as  one  travels  in  the  £'  direction. 

Figure  I  also  shows  the  coordinate  system  that  is  to  be  imposed  at 
each  point  of  the  prolate  spheroid.  The  prolate  spheroid  is  described 
by  the  equat i on : 


2.  2  2 

x  ty  +  =  , 

2,  2  ,  \  .22 

2  (u  -1 )  l  u 


which  generates  an  outward  normal  unit  vector  at  a  point  P  =  (x,y,z): 

2 

N  *  C  (x  ,  y ,  z  ^-y-)  (9) 

u 

where 

r  t  2o.  2.  2  /U^-1 » 2«-l /2 

C  =  (x  +y  +z  ( — — )  )  (1C 

u 

It  is  convenient  to  define  a  coordinate  system  at  P  where  N  defines  one 


axis,  and 


/  2  2 
x  +y 


-±  ■+  -+■ 

T  -  N  x  B 


define  the  other  two.  This  choice  of  £  is  stimulated  by  the  symmetry  of 
the  prolate  spheroid.  £  is  tangent  to  any  circle  defined  by  lidding  z 


constant  and  is  normal  to  N. 


,  jj m  um'iiiihtfiUii'Mh'-Sl'ii' 


6 


It  is  easily  seen  that  the  incident  and  reflected  rays,  k  and  k' 
respectively,  lie  in  a  plane  containing  N  at  P(x,y,z).  Without  loss  of 
generality  we  can  assume  that  both  It  and  it'  are  unit  vectors,  and  it'  can 
be  writ  ten 

it'  -  it  -  2  (it *N) N .  (13) 

Having  found  the  direction  of  the  reflected  energy,  it  is  now  neces¬ 
sary  to  derive  an  explicit  expression  for  the  spread  function.  Figure  2 
demonstrates  the  underlying  assumptions  regarding  the  spread  function. 

It  is  readily  seen  that  all  the  power  reflecting  from  the  rectangle  in 
the  BT  plane  described  by  the  sides  XT  and  y  B  must  pass  through  the  surface 
described  by  the  tips  of  the  vectors  t  k! .  Define 

R,  =  XT  +  tk,  *  -  titQ'  (14) 

and 


i$2  «  x!  +  tk2' 


(15) 


These  vectors  describe  the  sides  of  the  rectangle  at  the  tips  of  the 
tk' .  and  the  spread  function  can  now  be  defined: 


spread (tk' )  =  1 im 
U  Xh-0 
Y  -0 


V  •  VK2 

k  '  •  XTxyB 


(16) 


Equation  (16)  is  still  not  in  a  form  that  is  amenable  to  implemen¬ 
tation  on  a  computer,  but  it  can  be  made  so  very  easily  by  approximating 
k1,  by  the  expression 


dk  '  dk- '  ,  dkrt 1  dk  ' 

*  V  •  x  "d7~  *  dT~  ’  *  l}  *  1  (~dT~  *  — )] 

k  •  N 
O 

dV 

The  expre  ssions  i  nvol  ved  in  trying  to  eval  uate  — are  cumbersome 

— ► 

because  it  is  necessary  to  know  k'g  explicitly  in  terms  of  x,  y,  and  z, 


but  by  examining  Eq .  (13)  we  can  get 


■21*.$)*  -2(^)4!  , 


-2(k-f^)N  -  2  (k  -  N)  ^  . 
a  Y  d  Y 
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-¥ 

Equation  (9)  gives  N  in  terms  of  x,  y,  and  z,  and  -rr  and  —  can  be 

qa  a  y 

calculated  and  are  given  by  the  following: 


-*  2 
dN  r  T  T  u  - 1  x 
3—  »  C  ( T  ,  T  ,  T  — z— 
dX  x’  y  z  2 

u 

"c2[xVyTv+zTz(H^}  ]> 

y  u 


CCS  , B  , B  ii) 

d  Y  x  v  z  2 

1  u 

2  2 

-C2[xTx+yT  +zTz(ii-^i)  ]n. 
y  u 


Finally,  referring  to  Eq .  (6) 


%  ■ h  <vVr)  57  • 


assuming  that  D  =■>  t  and  that  L  (r) 
r  p 


1 ,  we  get 


spread (tk^1 ) 


It  will  be  noticed  that  Eq.  (21)  for  the  spread  function  involves  terms 
2 

in  I,  t  and  t  .  i  s  usually  examined  in  the  limit  as  t-*-».  0^  then 

becomes 


°b  =  1 


4r ! ’  ■  N  | 


t  dkQ'  dk^ 

ko  dx  *  dy 


3.0  RESULTS 


Aopendix  1  is  a  listing  of  the  program  that  was  used  to  find  the 

bistatic  cross-sectional  area.  SB  is  the  cross-sectional  area  as  given 

by  Eq .  (28)  while  SIGMAB  is  the  value  given  by  Eq.  (22),  retaining  the 
- 1  -2 

t  and  t  terms  in  the  final  expression. 

Appendix  2  is  a  graphical  representation  of  some  of  the  results  of 
this  p  rogram. 

Figure  3  shows  the  coordinate  system  used  to  orient  the  prolate 
spheroid  in  space.  The  major  axis  of  the  prolate  is  the  z  axis.  The 
transmitter  is  located  in  the  y,  z  plane  which,  because  of  the  symmetry 

involved,  does  not  affect  the  generality  of  the  program. 

Graphs  1-4  show  for  various  transmitter  positions.  The  observation 
plane  is  the  y,z  plane.  The  prolate  spheroid  has  major  axis  one  unit 
long,  and  minor  axis  .33  units  long.  is  in  square  units. 

Graph  5  shows  for  a  sphere  of  a  radius  one  unit. 

Graphs  6,  7,  8  show  o  in  varying  cases  when  the  observation  plane  is 

b 

"Ot  the  y,z  plane.  (These  are  all  for  a  prolate  spheroid  with  major  to 
minor  axis  ratio  of  3:1).  The  observation  plane  in  Graph  6  contains  the 
z  axis,  but  is  at  45°  with  respect  to  the  y  axis.  In  Graph  7,  the  obser- 
vat  ion  plane  is  the  x,z  plane.  In  Graph  8,  the  observation  plane  is  the  x, 
z  plane,  but  rather  than  having  the  transmitter  on  the  y  axis  as  in  Graphs 

7  and  8,  the  transmitter  is  at  30°  with  respect  to  the  z  axis  in  the 

y,z  plane. 

Graph  9  shows  when  the  major  axis  to  minor  axis  ratio  is  changed 
from  3:1  to  10:1.  In  all  other  respects,  this  is  the  same  as  Graph  2. 

The  observation  distance  was  equal  to  1000  units  in  all  the  graphs. 


Prolate  Spheroid 
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It  will  be  seen  immediately  that  for  a  prolate  spheroid  with  major 
axis,  1  unit,  the  maximum  value  of  is  tt.  Graphs  1-4  show  that  this 

maximum  value  in  the  y-z  plane  is  at  180°  minus  the  transmitter  angle  which 
is  exactly  as  expected  in  terms  of  reflections,  and  geometrical  optics. 

It  is  also  interesting  to  note  that  the  patterns  in  Graphs  1-4  are  exactly 
the  same  except  for  being  rotated  in  space. 

Graphs  4,  6,  and  7  show  another  aspect  of  this  pattern  in  space.  The 
pattern  has  maximum  width  in  the  plane  containing  the  transmitter  and  z  axis, 
and  becomes  increasingly  narrow  as  we  rotate  the  observation  piane  away 
from  the  y,z  plane.  It  is  to  be  noted  that  the  value  of  as  we  approach 
the  shadow  becomes  undeterminate.  In  the  x,y  plane,  *»  tt  for  all  angles, 
while  in  the  y,z  plane,  approaches  .02702. 

Graphs  2  and  8  show  the  same  behavior  as  Grap  .s  4,6  and  7,  except  that 
the  transmitter  is  now  at  an  angle  of  30®  with  respect  to  the  z  axis 
rather  than  being  on  the  y  axis.  Again,  the  pattern  approaches  being 
infinitely  thin,  and  is  indeterminate  as  we  approach  the  shadow. 

Lastly,  a  comparison  of  Graphs  2  and  9  shows  the  effect  of  increasing 
the  major  axis  to  minor  axis  ratio.  If  the  prolate  spheroid  becomes 
narrower  the  maximum  value  remains  tt,  but  the  pattern  becomes  much 
narrower. 

4.0  EVALUATION 

The  first  question  that  must  be  asked  is  whether  the  results  obtained 
have  any  meaning.  First,  the  monostatic  cross-sectional  area,  o  ,  for 
a  jrolate  spheroid  with  the  transmitter  on  the  z  axis  is  a  well-known 
quant  i  ty  j  and 
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a 

m 


ir(u3-))2*2 


2 

u 


(29) 


For  a  sphere,  a  *  it,  wh  i  ch  i  s  w  i  th  in  .  I  %  of  the  value  found  using 
m 

the  program.  For  a  prolate  spheroid  with  majo-  axis  to  minor  axis  ratio 

3:1,  (u  =  1.05,  2.  =  -95),  the  program  gives  o^  =  .02702,  while  Eq.  (29) 

gives  a  =  .02702.  In  fact,  for  ratios  3:1,  10:1,  and  22:1,  the  agree- 
m 

ment  was  exact. 

Second,  the  Crispin  bistatic  theorem  [8]  stipulates  that  given 

transmitter  and  receiver  directions,  k  and  k1  respectively,  in  the  limit 

of  infinite  distance,  a.  is  equal  to  c  as  seen  from  the  k  and  k‘ 

b  3  m 

direction.  Table  1  gives  a,  ana  a  for  various  transmitter  and  receiver 

b  m 

locations.  (The  angles  listed  are  with  respect  to  the  z  axis). 

As  can  be  seen  from  Table  I,  there  is  close  agreement  between  c, 

D 

and  a  .  The  small  discrepancies  shown  in  the  last  two  entries  are  easily 
m 

accounted  for  when  we  consider  that  since  the  distance  is  not  infinite, 

we  are  not  looking  at  the  same  point  on  the  prolate  spheroid  when  calcu- 

1  at i ng  o ,  and  c  . 

3  b  m 

We  can  conclude,  then,  that  the  program  is  giving  meaningful  numbers, 
and  that  these  numbers  do  represent  the  bistatic  cross-sectional  area 
of  a  prolate  spheroid. 


5.0  CONCLUSION 

The  prolate  spheroid  is  a  convenient  target  for  radar  studies 
because  it  is  convex.  F-om  a  geometrical  optics  point  of  view,  this  is 
convenient  because  it  implies  that  given  an  observation  point,  there  is  only 
one  point  on  the  prolate  spheroid  at  which  reflection  occurs.  As  such, 
there  is  no  need  to  account  for  the  phase  of  the  signal.  Any  surface  that 


1 4 

is  convex  at  each  point,  and  which  has  no  two  distinct  points  where  the 

outward  normal  vectors  are  in  the  same  direction,  can  be  handled  similarly. 

Hence,  the  program  written  could  very  quickly  be  changed  to  handle  a  figure 

that  is  ellipsoidal  in  all  the  major  planes,  or  even  a  parabola  of  revolution 

that  extends  to  infinity  in  one  direction.  It  is  simply  a  matter  of  being 

able  mathematically  to  describe  the  surface. 

When  the  outward  normal  vectors  are  not  distinct,  it  is  necessary 
i  k  r 

to  account  for  the  e  spatial  variations  in  both  the  incident  and 
reflected  fields.  This  can  be  done,  but  was  beyond  the  scope  of  this 
study  in  both  computer  funds  and  time  available. 

Such  a  program  would  be  a  major  step  forward  in  the  ability  to  find 
the  bistatic  cross-sectional  area  of  arbitrary  surfaces  in  space,  where 
it  is  not  clear  that  the  Crispin  Bistatic  Theorem  is  appropriate.  This 
would  be  an  area  for  further  study. 


TABLE  1 

A  COMPARISON  OF  o,  AND  o  FOR  THE  CRISPIN  BISTATIC  THEOREM 
b  m 


T  ransmi tter 
Loca  t i on 
(degrees) 


Rece i ver 

,  .  a, 

Location  b 

(degrees)  (units) 

_ _ (sq  uared) _  at 


0 

O 

-30* 

.02702 

60° 

1 

ON 

o 

o 

.02702 

30° 

150° 

3-126 

0° 

O' 

o 

o 

.04522 

30° 

90° 

.26481 

o 

m 

it  +  k‘ 
02  702 
02702 
I.  126 
.  04520 
.26464 


(The  observation  point  is  I0C0  units  away). 
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KEAL  L,«J,  THET  A,..*,PC3)  ,«i, E3TkP(3) 

REAL  K(3)»  6(3}»!U3)»T(3)fKPPlME(3)fpLf?(3) 

REAL  OGn ( 3) >  DLk P ( 3 ) , C 1 (3),C?C3),OGKP(3) 

REAL  FL(3),FG(3),SG,OIST,SIGMAd,OTSPNT(3) 

INTEGER  I, J#m, JP»IMCAY,Ii\iCA7 

LOGICAL  TEST,  RETEST 

REAL  MAG, ALPHA, BETA, KLIST 

PPIwT,  '  L  s  ' 

RE AD (5, 200)  L 
200  FORMAT (F 1 0 . 5 ) 

PR  If  T  ,  'Us' 

READ ( S  »  20  0 )  U 
PFINT,  ' Tm£ ( A  s  ' 

READ (5, 200)  TnETA 

KR IT £(03, 10a)  L,L, THETA 

lea  FORMAT  ( '  L  =  ',F6.A,'  u  s  THETA  r  ',F8.0) 

C  K  IS  THE  VECTOR  Thai  DESTINES  T^E  u I  RFC  T I  OfJ  OF  THE  INCOMING 
C  RADIATION.  Th£  AkGlE  Thf;T«  IS  fi  t  T  !.  C 0  AND  90,  INCLUSIVE. 

PI  =  ARCCS(-l.O) 

THETA  s  f’  I  *  1  HE  T  A  / 1  80 . 0 
K(l)so 

KC2)  s  «1.i)*S1h(1mKT.o) 

KC3)  s  « 1  . 0 ■* C O S  ( 1  H t  T  A  ) 

PRINT,'  DISTANCE  =  ' 

READ (5 , 201)  LIST 

20  1  FORMAT ('  ' , F 1 0  ,  )  ) 

PRINT ,  '  INCH  OF  A'nGLE  IN  V  =  (13)  # 

REA  u  (5, 202)  1/CAY 

PRINT,  '  I  iv  Ck  OF  A'-jGLL  I's  Z  =  ( 1  3 )  r 
READ  ( Si ,  20  2  )  IiCAZ 
202  FORMAT ( 13) 

RETESTs. FoLSt. 

DO  10  Jsl , 1  B  J , I:  CAT 
BETAS  J-l 

BETAS  PI*i»ETA/ibO.O 
DO  20  Ms  1,101,1 MC A i 
IF  (RETEST.  A‘iD.M.EN.  1 )  GO  TO  2^ 

IF('-.EU.l)  .<ETEST  s  .T'fuE. 

ALPHA  =  r-1 

ALPhh  s  P  1  *  ALH'h  u  /  1  ^  0  , 

C  CNN  V  R  T  TAKES  T  HE  OIsTa'CE  A -iu  Af  GLLF,  AND  CMP  VERTS  THESE  TO  HE 
C  PCIf  T  DP  SIRED. 

CALL  C'J.'.V'  r  (LIST  ,  ALPHA  ,:.!E1  A,CJ'-i,PNT) 
ft  PIT  E(03,  1  OB) 

F'V-ATl'  Tut  uHM-fVun.'*.  POT. IT  =  ' 


105 


TEST  s  .FALSE. 

c  if  the  puint  is  too  close  Tu  ihf  axis  that  k  is  om, 

C  THEN  THE  RULT1  'lES  9LO>v  DP. 

IF(  (MAGCObSHM  )-0!)T(0l;Sf',J.K))  .LE.l  .0E-1O)  TEST  =  .TRUE. 

.  IF( TEST)  GU  TO  40 

CsOBSPNTU  )**2  +  UflSPNT(2)**2+(*lW-l  .  A )  *  (Of*SPNT  (3 ) /U)  **2 
IF(C.IE.L**2)  TESTs. True. 

IF(  TEST)  aRITE(03, 107)  '  DBS  PUT  TOO  SMALL  ' 

.107  FORMAT  (  *  ',420) 

IF(TEST)  GQ  TU  20 

C  THE  PROGRAM  IMG*  GUES  AGUUT  l.vPUSr'G  AJ  N,P,T  FRA E  ON  The 
C  SPHEROID  AT  THE  PUluT  OF  REFLECTION.  INC  IDF'i  T  At  LY ,  IT  MUST  FIRST 
C  FIND  THE  POINT  OF  REFLECT  I  Of!.  *EMCE,  FINDKP. 

CALL  NURH(OdSPNT,Cl) 

CALL  suncci,K,ci) 

CALL  NURfi  (C  1  ,  -i) 


CsL*(U**2-l«0)/SGRT((o**2-1.0)*(,<(l)**2  +  N(2)*‘,2)«‘U**?*M(3)**2) 
P  C 1 )  s  C  *  N  C  1 ) 

PC2)sC*N(P) 

P(3)sc*U**2/l’J**^-l,0)*  :  (3) 

CALL  Fli  UKPd-  , N,D,T  ,KPiiIf-L) 

CALL  SUIUuBoP  iT,P,ESTSP) 

CALL  NO R rt ( E 3  T  K P , E  S  T  k  P  ) 

CALL  ADOCKP-vi  '.E,  E3TKP,C1  ) 

CALL  nUk'UC  l ,  kPkI  -ie) 

CALL  SUP(rPHI'1E,K,rj  ) 

CALL  NURMCl,  O 

CALL  SUPCnfRI  LSTn  P  ,  f  T  m'  ) 

IF  (i-'.AG(ESTKP)  ,GT.  i  ,0E-'4)  GO  TO  1  Oh 


CALL 

CALL 

CALL 

ir<i-U 


IF(LOT((  ,-O.GE.0.0)  TESI  =  .T''UF. 

IF  (TEST)  v-.HlTt  Co3,  1  0  7)  '  1  o  T  !,C  SHAl'n 
If  (1 EST)  GO  TU  2  0 
CALL  SOU ( DtiSH JT , E , C J) 

K  D  I S  T  =  O  A  G  (  C  1  ) 


IT  IS  NECESSARY  TU  ucl  t'O-.  t  FE  A  TURPS  OF  THE  CURVATURE  OF  THE 
SPHEROID  AT  T>.l  PUI  .1  OF  i.'EfLCn  ION.  T HE  DERIVATIVES  OF  N  ARE 
TAKEN  I"  T  up  R  A  i,i)  T  ■)  I  !'EC  1  I  U'<3 .  LA".'’<f>A,  A-,0  G  AM**  A  APE 
DLMNY  COORDINATE  uESCr  1  ?  |. >i  $  1  M  T-E  T  \Ci  OJRFCT  IONS. 

CALL  D  E  h  i  v  (  P  ,  ,  .i ,  L  '■ ) 

CALI  DE^I v ( P, T , u»  OG  .  J 

If  (GUT  (  0  L :  ,  -')  .GT.u.O)  CALL  SC  L  :L  T  f-1 . 0 , 01.  >>,  f  ■  L  ’ :  ) 

1F(I  OT  C  iU  ,  1 )  . G  I  .),!»)  C  ALI  SCL  'L  T  (-1  . o  ,  nr,rj,  crp.. ) 


IMlikhUililtUll 


C  USING  ULN  A \ 0  UGn,  IT  IS  PUSSlaU-  TO  CALCULATE  THE  DERIVATIVE  OF' 
C  KPHIWE,  Th£  VECTOR  FK<»‘-  T HE  PUI*  T  '.in  THE  P.S.  TO  THE  DHSFMVATION 
C  POINT,  IN  THE  ThO  D I  k£C  T I  M„S ,  AND  THEN  TO  FIND  THE  SPREAD  OF  THE 
C  RA0IAT1ON  AS  IT  TRAvEt  3  AaAY  FK-Im  tl't  P.S. 

CALL  SCLMLTCL'UT  C*,*  )  ,DLf  ,C1) 

CALL  SCLMLTCOUT  (OL'J  ,<),,;,  C2 ) 

CALL  ADC(Cl,C2,DLKP) 

CALL  SCLNLT 1-2.0, DLL P,DLKP) 

CALL  SCL-LHUOT(w,M  ,r.6\,Cl ) 

CALL  SCL  «LT(  OUT  (DG.l,*)  ,-n,C2) 

CALL  A0DCC1  ,C2,L'GsP) 

CALL  ADD (C 1 , C2, DGrP) 

CALL  SCLMLT(-2.0,DuKR,0GKP) 

CALL  CROSSlo,i)GSP,Cl) 

CALL  CROSSIDLnP, T,C2) 

CALL  AL>0CC1,C2,CI) 

CALL  SCL-'l  T  (  I.O/LDIST  ,C1  ,C1) 

CALL  SCL!,L  T  ( 1  ,0/VL}lsr**2*'.,C2) 

CALL  Al)D  (C  1 ,  C2 ,  C  1 ) 

CALL  CROSSCOLaP, UGKH/C-?) 

C  SB  IS  THE  HISTATIC  C'OJSSSEC  TlOnAL  A  WE  A  IF  wf  IGNORE  THE  D  AND 
C  0**2  TER  *'5  I  N  Th£  SPR'Em'  EXRr  F.SSl'W',. 

Sfis4.o*PI*UUT  (  h  ,  h  H  0 1  -if  )/D.i  f  ( C  2  ,  K  P 55 1  E  ) 

IF(StT.LT.O.O)  S  ■»  =  -l.0*3:< 

CALL  AL»D(C1,C2,C2) 

C  SIG^'AF*  IS  The  True  CLUSSSECTJOnAL  AREA, 

SIC'-Ats-l .  0*21  *0  IT  (.*,  hPHI  'E  )  />: 01  (  C2,  KPRI-’E) 

IF  (S  ICHAH.LT.  0 .0)  S  I G'  A  <>s-  1  .  0  *S  I  G  •’  A  n 
ARirt  (03,102)  'Sh  =  ',bn,'  Sir,.!ip  =  * ,  S  IG!'AH 
102  F  C'  R  A  T  C '  ',  A5,F  10.5,  -U  1  ,F  J  O.S) 

20  CONTINUE 

10  CONTINUE 

STOP 
E  f-1  D 

HEAL  FlP.CTM'v  vagi/.  ) 

REAL  A (3) 

M  AG  =  Si»kT(A(1)**2  + ».(<>)*♦  £♦  A  (3)**?) 

RLTuKii. 

END 

SUHROUT  INF.  i  U  R  C  A  ,  r_  ) 

HEAL  A(3),RC3),  C,  'AF- 
CsSL’RT(A()  )*  *2+.»  (,*)+*.?♦  <i  (3)**?) 

IFCC.LE.l,  i»E-cO )  .-.KITE  (03,12)  '  DIVISMw  DV  ZERO  ' 

12  FORVATC'  "  , n 1  -  ) 

C=1 . 0/C 

CALL  SELL  I  l  C  ,  A  ,  i  • ) 

H  t  T  L  H  N 

ELD 


S  JBHOuTINE  FI*i)ivFlX,N,HiT'kPf'TMn 
REAL  K(3)»N(3)icT(3),T(3),‘\pcIMEf3) 
REAL  Cl(3),C2(3),MAb 
b  ( 1 )  s  ( N  ( 2  )  ) 
b(2)*-w(l ) 

B(3)sO.Q 

IF(f  AG(^) .LE. 1 . 0h-20 }  B  (  1  )  =  1  ,  0 
CALL  »NOR:'UB/b) 

CALL  C«OSSU,U,  T) 

CALL  SCLMLT(-2.0*u0T(K»-)#'J#CI) 

call  AU'>(K,cif KPbr-e) 

RETURN 

EM) 

SUBROUTINE  ADOCa,m,C) 

REAL  a(3),  d(3),  C ( 3 ) 

INTEGER  I 

l)U  10  1  =  1,3 

c(n  =  A(i)t8cn 

CONTINUE 

RETURN 

EM) 

SUBROUTINE  Sud(AfOfC) 

REAL  A(3),B(3),C(3) 

INTEGER  I 
DO  10  1=1,3 
C (I) =  A ( I ) - b  C I  ) 

CONTINUE 

RETURN 

END 

SUBROUTINE  SCL'-'LT  (A,B,C) 

REAL  A,i'C3)  ,CC3) 

integer  1 

UG  10  1=1,3 
C(I)aA*il(I) 

CONTINUE 

RETURN) 

END 


REAL  FUNCTIJN  iO)T(A,N) 
REAL  A(3),tH3) 

INTEGER  I 
D  U  T  =  0 

00  10  1=1,3 
cm  =  DOT  ♦  » (  I  )*<()) 
CONTINUE 
m  E  T  OR 


SliywUU  T  i  ''■'t  CRUSS( A,P,r) 

RE AL  A(3).U(3)»  Cl  3) 

C<l)sA(2)«d(3)-A(3)*»H2> 

C(2)»A(3)*«(l)-A(l)*R(3) 

C(3)sA(l)'*hl2)-AC^)*l.  Cl  ) 

RETURN 

END 

SUBROUTINE  CUU'rtl  (A,B,C,0) 

REAL  A,A#CiI>(3) 

REAL  P,fv 

LOGICAL  test 
TEST  s  .FAL3E. 

D(1 )=0.C 

0(2)  s  A* SI •'»(*) 

0(3)  =  A*CUSCB) 

IF(SIN(C)**2.LT.1  .VE-1  0.-Vji>.C  »S(C)  .01  .f'.'^c-eriiPi 

IF(SIN.(C>  W.LT.l.OE-W)  »EST  =  ,T«UE# 

IK (TEST)  0(Da-A*SI\(,-.) 

IK  (TEST).  ,0(3)8-a*C;1S(-) 
imtST)  RtTl-Ri. 

IF  (Sir  (  .  )  *  *  2  ,  L  T  ,  l  ,  A  E  -  l  n  )  RET1 1  c  . 
1K(C0S(,’)**2.LT.1  .  E  •  l  •'■ )  p  =  r.,o 
IK(C0S(C)**2.LT.l,0f.-u>)  =  A.n 

IF(Cl)S(B)**2.GT.l  .oe-u;)  P  s  i.*/ian(».o 
IF(CUS(C)**2.GT,  1  ,i*C-tO)  0  =  1  .  0/T  A-.(C) 
0(1)=A/S'JR1  (1.0t(.J)**2t-(P)W) 

0(2)80(1) 

0(3)80(1) *P 

RETURN 

END 

SUrtKOUTI-JE  UEKlvtPrl  (UP’L'O 
REAL  P(3)  t  T  (  3  )  »U,i)L-.C3) 

REAL  A  (3)  pU3)  #C  (3)  ,(>(  5)  #  -ap 
A(l)sPU) 

A(2)8P(2) 

A(3)=P  (3)  *  (<J*  *2-1  .  0  )  /  :*  *2 
B(l  )sT(D 
6(2)sT (2) 

B(3)aT  (3)*('J**2-l.'»)/ii**2 
CALL  SCLNL1  (1.0/*.AG(A),M,C) 

CALL  SCL"M.T  (OUT  (/»,*),  A,  n) 

CALL  SCL'VL  T  (1  ,')/■■  hi*  ('•)**  3  ,  »  ,  0  > 

CALL  S'JR(C,i.’,»L  ■) 

R  E 1  UR- 


j t»r"'  Mpm  rut  mmm  «t«w  *v*w 


inr^1  fnpif  < 


Y  Axis 


Graph  2:  Cb  for  a  Prolate  Spheroid  with  the  Transmitter  at  30°w.r.f . 
the  Z  Axis  (major  to  minor  axis  ratio  3:1) 


Y  Axis 


Graph  3:  Ofa  for  a  Prolate  Spheroid  with  the  Transmitter  at  60°  w.r 
the  Z  Axis  (major  to  minor  axis  ratio  3:1) 


Z  Axis 


aph  6: 


for  a  Prolate  Spheroid  with  the  Observation  Plane  at  45° 
(major  to  minor  axis  ratio  3:1) 


X  Axis 


Graph  8:  for  a  Prolate  Spheroid  with  the  Transmitter  at  30°  w.r.t. 

the  Z  Axis,  with  the  Observation  Plane  being  the  X,  Z  Plane 
(major  to  minor  axis  ratio  3:1) 


Z  Axis 


Graph  9:  for  a  Prolate  S 

with  the  Transmit 
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